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Abstract 



Radiative corrections due to pair production in process of electron-positron annihi- 
I lation are considered. The main attention is paid to the corrections of the third order. 

Simultaneous emission of photons and pairs is taken into account. The leading log 
O . approximation and convolution procedure are used. Exponentiation is discussed. Theo- 

in \ retical uncertainty in description of higher order secondary pairs is estimated. 

' PACS: 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 

i!B ■ 1 Introduction 

Qh! The high statistics and accuracy of experimental measurements at LEP and SLC result in the 
outstanding level of precision in definition of the Standard Model parameters (see Ref. [H |2] 
and references therein) . In this situation effects of radiative corrections are of great importance. 

One of particular contributions to the process of electron-positron annihilation is the ra- 
diation of secondary pairs. In comparison with the photon radiation, it is relatively small, 
because it appears only starting from the 0{a'^) order. Nevertheless, the total effect of pair 
production can reach several permil (or even e few percent in particular cases) and should be 
taken into account. We will consider here the reactions of the following type: 

e+ + e- ^ / + / + (/'/'(7)), (1) 

where // is the so-called primary fermionic pair, /'/' is the secondary one, and (7) stands 
for the usual photon radiation, which can take place in our case starting from the third order 
O (a^). The primary pair is the one registered in the particular measurement. The secondary 
one can be either real or virtual, and has to be accounted as a radiative correction according 
to certain experimental conditions (with some restrictions (cuts), if required). In principle, 
the pairs can be of the same flavour, like in the process e~^e~ —>■ 2fi~^2fi~ . When both pairs 
are seen in the detector, one has to apply an algorithm to choose the primary pair. 
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The complete second order QED calculation for production of the initial state secondary 
pairs (e+e~, fi^fi~, and hadrons) was performed in Refs. jnill]. The corresponding final state 
effect was considered in Ref. jS]. 

The interference between amplitudes for the initial and final state pair production has not 
been calculated systematically. An estimate of this effect will be presented in the Conclusions. 

It is observed, that the 0{a^) approximation is not adequate to provide the desirable 
precision, especially in the region close to the Z-peak. From another side, at LEP2, when 
the radiative return is allowed, QED radiative corrections reach several hundred percent and 
should be evaluated with care. In order to achieve the precision, required by the experiments 
(see Ref. one should consider higher orders of the perturbation theory. Another important 
point is that in the Z-peak region the sub-leading corrections are comparable in size with the 
leading ones. So, one has to go beyond the leading log approximation (LLA). 

The goal of the present paper is to describe and discuss the treatment of pair corrections, 
implemented in the semi-analytical code ZFITTER v . 6 . 30 |6i| . The main attention will be paid 
to the third order corrections. Some numerical illustrations will be given. The uncertainty of 
the approach will be analyzed in the Conclusions. 

2 Pair corrections in higher orders 

Starting from the 0{a^) order, one has amplitudes with simultaneous emission of pairs and 
photons. We treat them as a part of pair corrections, because they have not been taken into 
account in the traditional photonic RC. The main subject of this article is to analyze this 
effect. 

As usually in e~^e~ annihilation, the most important contributions are from the initial 
state (IS) radiation (real and virtual). In a typical high-energy experiment, the conditions of 
the Kinoshita-Lee-Nauenberg theorem are fulfilled for the final state (FS) radiation. That is 
why the large logarithms (mass singularities) cancel out in the FS correction, and the latter 
becomes relatively small. From the other hand, a specific energy dependence of the kernel 
cross section can provoke an enhancement of the IS correction, like it happens at the Z-peak 
or in the radiative return to the resonance at LEP2. 

Pair corrections to the process of e^e~ ff become ambiguous, when one observes 4- 
fermion final states. Nevertheless, the bulk of the pair RC to the 2-fermion observables is 
coming from the kinematical situations, when the secondary pair is either virtual or escapes 
registration (in the beam pipe or because of a small energy). On could try to subtract from 
the experimental data all real secondary pairs by means of a Monte Carlo generator. But that 
would destroy the cancellation between the real and virtual pairs, and the corrections would 
become even more pronounced and dangerous. One can get a serious numerical instability, 
trying to subtract both the virtual and the real pairs by means of different codes. Even if the 
subtraction of the pure second order pairs is evaluated perfectly, we have to keep in mind the 
higher order corrections due to instantaneous emission of pairs and photons. The cancellation 
between the virtual and real pair RC is important not only in the O (a^), but in higher orders 
as well. Indeed, we would end up with a completely wrong third order corrections (the error 
would exceed the entire effect), if we were using for the subtraction of virtual and real pair 
contributions two different codes, one with the third order pair-photon effect and another one 
without it. 
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Figure 1: Different types of Feynman diagrams for real pair production. 

Here we will follow the definition of signal and background according to Feynman diagrams, 
as discussed in Ref. In Figure 1 we show several representatives of Feynman graphs for 
real pair production processes. Only the diagrams Fig. la and Fig. lb (plus the ones with pair 
production via a virtual photon emitted from the initial state positron of the final state anti- 
fermion /) are to be accounted as the signal. These diagrams describe the non-singlet (NS) 
mechanism of pair production. The diagrams of the type (a) and (6) with the production of a 
secondary pair via a virtual Z-boson are treated as background. The so-called multi-peripheral 
process is described by Fig. Ic. Diagrams of type (d) represent the singlet mechanism of pair 
production. Their contribution is considered as a background too. The pure electroweak four- 
fermion processes with ly-boson intermediate states fall into a separate class from both 
the theoretical and the experimental points of view. The relevant (signal) Feynman graphs 
for virtual pair corrections are given in Figure 2. 

In other words we will count only the non-singlet secondary pairs produced via virtual 
photons Fig. la, 6 (some numbers for the singlet pairs will be presented ). The contribution 
from the so-called background diagrams, like multi-peripheral, double-resonant (ZZ), and 
singlet pair production, are supposed to be subtracted from the data. Interferences between 
the signal and background amplitudes are typically small, but some of them, like the ones 
between the IS singlet and non-singlet pairs, should be taken into account. We will not use 
any direct cut on the phase space volume of real secondary pairs (an inclusive treatment of 
radiative events). Only the condition (see Eq. Q) on the primary pair invariant mass s' will 
be applied (the dependence on the secondary pair invariant mass will be shown for the FS 
correction) . 
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Figure 2: Feynman diagrams for virtual pair corrections. 
2.1 The final state pair RC 

In the lowest order for the FS pair corrections we use the results of Ref . |3| . Our aim here is to 
estimate the effect of higher orders. Concerning the interplay with the IS photonic corrections, 
it is claimed, that the final state pair corrections should be realized in the multiplicative way: 

O- = aBorn(l + (5lS7)(l + (^FSP), (2) 

where Sis-y stands for the initial state photonic correction, and 5fsp gives the final state pair 
one. At LEP2 energies, when the radiative return to the Z-peak is allowed, we have huge 
values of Sis-y (typically above 100 percent), and the additive summation of the corrections 
(1+5is^ + 5fsp) fails completely. In principle, the initial state pair correction should participate 
in formula (j2)) together with the photonic one, but, as one can see from numerical illustrations, 
an additive treatment of the IS and FS pair corrections is acceptable as well. 

The effect of simultaneous photon and pair radiation from the final state can be estimated 
by multiplying the corresponding contributions: O (a) FS photon RC and the O (a^) FS pair 
RC. One can see that this source does not give more than 0.05 permil in the accepted signal 
definition. This estimate was confirmed also by a computation (analogous to Eq. (jS))) of the 
O [a^) LLA result for the FS non-singlet pairs. It is worth to note here, that additional 
strong cuts on the final state photons and pairs can increase the FS corrections and revive the 
problem of the corresponding missing higher order terms. 

In Tabled we give the value of correction (5fsp due to the final state pair production in the 
process of annihilation into muons as a function on cuts on the invariant mass of the primary 
and secondary pairs; 

-^primary = ■S > ^min-S, -^secondary — '5-Pcut, 

W = 10'(2|?^-l]. (3) 

where crJsp ^ is the cross section with the photonic and IS pair corrections taken into ac- 
count. One can see that the dependence on the invariant mass of the secondary pair is very 
weak. That confirms our expectations, because the contribution of large-invariant-mass pairs 
(MggPQj^^g^j.y ^ m'j, where is the mass of a pair constituent; ruf = for the case of hadronic 
pairs) does not contain any large logarithm, and so it is suppressed by the relative factor 
The same should be valid for the IS pairs: the relative contribution of large-invariant-mass 
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Table 1: Radiative correction (5pgp ^ in permil for different cuts. 

secondary pairs is negligible. This statement concerns, of course, pair production via virtual 
photons, while the double-resonant ZZ-mechanism (or VTiy) has to be treated separately. 

2.2 The initial state pair RC 

The combination of the initial state pair correction with the final state photonic one should 
be treated in the same way as Eq Q, while it is not so important numerically. We will 
concentrate now on the main problem, when both the pair and the photons are produced in 
the initial state. 

The initial state pair contributions to the corrected cross section is presented as the con- 
volution of the boosted Born cross section with the so-called pair radiator: 

(5crP^'''= r dzd{zs)H{z) = a{s)H^+ T ^ dz a{zs)He{z). (4) 

Where ifA represents the impact of virtual and soft (small-energy) pairs, A is an auxiliary 
parameter {ml/ s <C A -C 1); numerical results should be independent of its value. In higher 
orders, when we have instantaneous emission of pairs and photons, we will still call such 



5 



corrections as the pair ones. In our computations the O (a^) formula for the pair radiator is 
taken from Refs. [Sll^- 

Going beyond the lowest order, one starts usually with the leading logarithmic approxima- 
tion. Indeed, the main contribution to the radiative correction under consideration is coming 
from the electron pair production, because it is reinforced by powers of the large logarithm 
Le = ln(s/mg) (the corresponding leading logs for pairs of other flavours can be constructed as 
well). There are two different mechanisms for electron-positron pair creation in our process: 
the singlet and non-singlet ones. The first one is absent for non-electron pairs. 

At the Z-peak it is worth to try to re-sum higher orders, because the convergence of 
the perturbation theory there is slow. The first exponentiated formula for the IS LLA pair 
production was suggested in Ref. [7j. The process of one pair production was supphed there by 
emission of an arbitrary number of photons. This approach gives a reasonable approximation 
for the leading corrections close to the Z-peak. But it does not include some significant next- 
to-leading terms, and even the known third order leading logs are not reproduced completely. 

Another way fH] is to use the Mellin transformation to solve the QED evolution equation 
for soft and collinear photonic contributions supplied by the running QED coupling constant. 
The authors received a simultaneous exponentiation of photonic and pair corrections, which 
takes into account important next-to-leading terms. Unfortunately, their formulae require 
additional evaluations (inverse Mellin transformations). Moreover, the applicability of their 
method to the mixed pair-photon corrections is not clear. 

In Ref. |9j a phenomenological formula for simultaneous exponentiation of photonic and 
pair radiation was proposed. The correspondence of the exponentiated formula to the pertur- 
bative results was shown for the case of real hard radiation. Nevertheless, the structure of the 
radiator function, suggested in [9^ , does not allow to verify analytically the correspondence for 
the soft and virtual part of the corrections. 

2.2.1 The leading logarithmic approximation 

Let us start our evaluations with the leading logarithmic approximation. It will serve us as 
a guide. The third order LLA radiator function can be taken from the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi evolution equation (we dropped below the pure photonic corrections): 



^^(3)pair ^1 dza{zs)(^^{Le-l 

where 

R^{z) = 
is the singlet splitting function; 
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are the non-singlet splitting functions. The convolution of the singlet splitting function i?^ 
with the first order non-singlet one is known [TUj : 



+ 2 ln(l - z) ]R'{z) + (1 + z){- In" z + 4Li2(l - z)) 



+ -{-9-3z + 8z^)lnz + -[ 8 + 8z + 3z 



(9) 



2.2.2 Beyond the LLA 

But, as we learned from the O (a^) level from the complete calculations, the contribution 
of the sub-leading correction in the region close to the Z-peak is very large (say, about one 
half of the leading one). In order to account the most significant part of the sub-leading 
corrections we consider the convolution of the 0(a;^) pair radiator with the ordinary 0{a) 
photonic radiator, proportional to the P^^^ splitting function. In this way we receive the main 
part of the 0{a^) leading logs, proportional to P^'^\ and the sub-leading terms enhanced by 
\n{l- z)/{l- z), like \n{l - z) / {1 - z) and L\n^{l - z) / (1 - z). Note that the convolution 
as well as exponentiation can not give the true complete sub-leading formula. And, contrary 
to the case of pure photonic radiation, they are unable even to cover the full leading log result: 
the term with P*^^^ is not reproduced (this term however has a small coefficient and no any 
ln(l — z), so it is negligible numerically in our task). This term corresponds to the physical 
situation (Feynman diagrams), when the extra photon is attached to the secondary pair. 

In fact the convolution gives a part of sub-leading terms coming from the kinematical 
regions, where both the pair and the photon are either soft or emitted collinearly to the 
parent particle momentum, while there are other sources for the corrections, like, for instance, 
emission of a collinear pair and a large-angle photon. But we believe, that the main terms with 
reinforcements are reproduced correctly: that follows from the general experience in leading 
log calculations. Moreover, for the case of pure photonic radiation this has been checked by 
direct perturbative calculations. 

In this way we received the expressions for the third order pair radiator. For leptonic 
non-singlet pairs we have 
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The large logarithm Lf = \n{s/m'j) stands for different lepton flavours: f = e, fi, r. The 
corresponding soft+virtual contribution is given by 
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The cancellation of the auxiliary parameter A was checked analytically and achieved numeri- 
cally. The largest contribution appears at the Z-peak, where the integral over z is practically 
limited by the Z-boson width: (2;mm, 1) ~^ (1 ~ ^z/Mz,l)- And therefore integrals from 
terms, like ln(l — z)/{l — z), are large. We checked that for the real hard emission there is 
an agreement between the most important terms in Hq {z) and the corresponding terms in 
expansion of the exponentiated formula from Ref. [HI. Such a correspondence between the 
exponentiation and convolution procedures is well known also in the case of pure photonic 
radiation. 

A similar expression was derived for the case of hadronic pairs: 
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(13) 



Quantities -R(oo), -Ro,i are defined in Ref. In contrast to Ref. [0] we extended the hadronic 
pair contribution to the third order dynamically, but not by a static coefficient. That is really 
important at LEP2 energies. 

The LLA leptonic contributions, which were not reproduced by the convolution (the terms 
with P(i) and P(^) (g) R' functions in Eq. ®), should be added to Eqs. (UHl), (HH). The third 
order effect appeared to be not small numerically. So, we have to look also at the fourth order. 
The corresponding contribution is estimated by means of the leading logs ^T] (non-singlet 
channel only): 



12 ^ ^ 216 ^ ^ 108 ^ ' 



(14) 



Only the term with the P*-^-* splitting function ^Oj seems to be numerically significant. It 
agrees with the corresponding term in the expansion of exponentiated formula in Ref. [HI . The 
size of the fourth order, which is estimated from the above formula for electron pairs, is rather 
small. So, we can be less rigorous: in order to estimate the fourth order contribution from 
other kinds of pairs and the large sub-leading terms, we drop them simultaneously. That was 
checked in lower orders to be a reasonable approximation. So, in the 0{a^) we keep only the 
leading logarithmic formula (fT^ for non-singlet electron pairs. 

In Table El we present the results for the different contributions of pair corrections to the 
process of annihilation into hadrons. The value of correction due to pairs is defined in respect 
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+ 16.44 



Table 2: Different contributions to 5fgp in permil. 



to the cross section for annihilation into hadrons with pure photonic corrections taken into 
account: 

fcp = 10^(^-l)- (15) 

The cut-off on the primary pair mass s' > Zmin ■ s, 2;min = 0.01, acts in the same way on real 
photons and secondary pairs. 

In Table El we present the values of 5jgp ^ for different cuts and beam energies. In 
the first two columns we give the contributions of singlet pairs and their interference with the 
non-singlet ones. The latter, in principle, should be counted as a part of the signal. The 
third column contains the NS contribution (our best estimate with all pair flavours summed 
up). It can be compared with our implementation of exponentiated pair corrections. Note, 
that in both the last two columns the contribution of hadronic pairs is not exponentiated, 
but it is supplied according to our approach (we used result ^ supplied with our third order 
contribution (IT^. 

For the column "|7j", the known complete second order result was supphed by the NS 
higher order part of the exponentiated formula, derived in this paper. The discrepancy with 
our best estimate originates from an improper reproduction of the third order correction by 
the exponentiation in Ref. j7]. Note that the largest difference appears at LEP2 energies 
for small 2;min values, which means the dominance of the effect of radiative return to the Z- 
peak due to hard pair emission. While the exponentiation [7j is well justifled only for soft 
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radiation, which was the only one relevant at LEPl. The last column was calculated by 
using the exponentiated formula of Ref. [H]. The approximation for the higher order hadronic 
contribution (static coefficient before the second order contribution), suggested in that paper 
has sense only at LEPl energy on the top of the peak and the left wing; everywhere else it 
fails. Errors in the third column represent our estimate of theoretical uncertainties (see the 
Conclusions) . 

As could be seen from Table El the contribution of singlet pair production becomes signif- 
icant for small values of z^^i^. In data analysis such events are supposed |2] to be extracted 
from the data together with the two-photon (multi-peripheral) process [T^. We underline, 
that the procedure should be accurate and well understood, because the events with singlet 
pairs and multi-peripheral production have completely different kinematics and signatures in 
the detector. At LEP2 energies the contribution of singlet pairs becomes really important, if 
the returning to the Z-peak is allowed (for ^ 0.25). The large contribution is coming 
from the peculiar 1/z term in the function R^, Eq. Note, that this function is multiplied 
by the kernel cross section which contains z~^ too. In this way the second order correction due 
to collinear emission of singlet pairs can give several percent. This did not happen at LEPl, 
because the kernel cross section outside the Z-peak is relatively small. 

The difference between the column "non-singlet" and the last one "P]*" is coming just 
from the difference of the perturbative and exponentiated treatment of higher order leptonic 
pair corrections. The difference is less on the right wing of the Z-peak, where the real radiation 
is dominant, because it provides returning to the peak. For real radiation we expanded the 
exponentiated formula and saw an agreement in the main terms. On the other side of the 
peak (and on the top) the virtual and soft corrections are dominant, because any real hard 
radiation leads to a huge reduction of the kernel cross section in Eq. (j3]). In this case we have 
to compare the soft+virtual part of our radiator Ha with the integral 



where p{z) is the exponentiated formula P, and p^{z) is the pure photonic exponentiated 
radiator ^3] • The latter is in perfect agreement with the perturbative results both in the hard 
and soft+virtual regions, while the formula with simultaneous exponentiations of photonic 
and pair radiation has a very peculiar structure for z 1, and so Eq. (fTB)) can be hardly 
expanded in a series in a to compare with the perturbative results. The problem is seen in 
the region of very soft photons {z —>■ 1) below the threshold of pair production, where the 
part of the exponentiated formula, which describes pairs, still represents a sharp dependence 
on 2;-value, instead of giving a constant coefficient from virtual pairs. 

It is worth to note the weak dependence on the Zmm parameter of the corrections at LEPl 
energies. So, the bulk of the corrections is given by soft and virtual pairs, which are not mixed 
during the observation with the primary pair from electron positron annihilation. 

3 Conclusions 

First we have to remind about the effect of the initial-final interference (IFI) in pair corrections. 
We suggest to estimate it by using the known effect for the IFI in photonic corrections. The 




(16) 



l-A 
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£^cm [GeV] 


singlet 


SxNS 


non-singlet 









Zrmn = 0.01 


88.1867 


+5.61 


-0.06 


-1.96 ±0.28 


-1.52 


-1.76 


91.1867 


+0.81 


-0.01 


-2.48 ±0.22 


-2.16 


-2.31 


94.1867 


+2.78 


-0.07 


+0.22 ±0.41 


+0.31 


+0.16 


189 


+51.03 


-1.08 


+ 11.78 ±0.82 


+ 14.17 


+11.39 


206 


+56.69 


-1.13 


+12.17 ±0.88 


+14.77 


+11.74 




-^min 0.1 


88.1867 


+0.24 


-0.03 


-2.17 ±0.28 


-1.78 


-1.94 


91.1867 


+0.03 


-0.01 


-2.51 ±0.22 


-2.20 


-2.33 


94.1867 


+0.12 


-0.05 


+0.12 ±0.41 


+0.19 


+0.07 


189 


+12.48 


-0.91 


+10.80 ±0.81 


+ 12.81 


+10.53 


206 


+18.14 


-0.95 


+11.13 ±0.86 


+ 13.40 


+10.83 




2^min 0.5 


88.1867 


+0.00 


-0.01 


-2.41 ±0.28 


-2.04 


-2.17 


91.1867 


+0.00 


-0.01 


-2.54 ±0.22 


-2.23 


-2.36 


94.1867 


+0.00 


-0.04 


-0.00 ±0.41 


+0.06 


-0.05 


189 


+0.02 


-0.05 


-1.05 ±0.09 


-0.80 


-0.96 


206 


+0.02 


-0.05 


-1.09 ±0.08 


-0.83 


-1.00 




•^min 0.9 


88.1867 


+0.00 


-0.01 


-2.58 ±0.29 


-2.20 


-2.33 


91.1867 


+0.00 


-0.01 


-2.59 ±0.22 


-2.28 


-2.41 


94.1867 


+0.00 


-0.04 


-0.32 ±0.41 


-0.25 


-0.37 


189 


+0.00 


-0.02 


-2.33 ±0.15 


-2.03 


-2.22 


206 


+0.00 


-0.02 


-2.40 ±0.16 


-2.10 


-2.29 



Table 3: Radiative correction 6 for different cuts. 

bulk of the IFI pair correction is provided by the terms of the order O (a^L) (the sub-leading 
logarithms). The terms are coming from the kinematical region of small-invariant-mass real 
pairs (and from the corresponding region of integration over virtual pairs). Thus, one can 
use the factorization of the correction factors, corresponding to the photon emission and the 
pair production. So, we just multiply the IFI photonic correction by the factor aL/(37r) and 
receive the estimate for the IFI pair RC [L = ln{s/ml ^^)). The typical value of the photonic 
IFI RC at LEP energies is a few permil, therefore one can safely neglect the corresponding 
pair contribution. Note, that in this way we considered only the non-singlet mechanism of 
pair production. The interference of the initial state singlet pair production amplitudes with 
the final state NS ones is expected to be of the same magnitude. 

Let us discuss now the uncertainty in the pair corrections. We mean just the indefinite- 
ness in the treatment of hadronic pairs and of the O (a^) (and higher order) leptonic pair 
corrections. Other sources, relevant for concrete experiments, like the problem to distinguish 
between the primary and secondary pair or the proper subtraction of the background, have to 
be discussed elsewhere j2]. 

First we suppose, that 15% is a reasonably safe estimate for the error in the second order 
hadronic pair contribution itself. In this case we use a phenomenological approach based 
on the measurements of electron-positron annihilation into hadrons at low energies. The 



12 



absence on the market of a plug-and-play parameterization of the i?-function (ratio of the 
annihilation cross sections into hadrons and muons) does not allow to put a better precision 
level. For our computations of the FS hadronic pairs a simple parameterization of the function 
was constructed. The dispersion relations were used to look at the corresponding vacuum 
polarization effect and check the function. The stability in respect of the fluctuations in 
involved data was verified. A good agreement was observed with the published numbers on 
hadronic pair corrections 01 E] in both the initial and the final states. 

The uncertainty of the higher order leptonic pairs is derived as follows. We assume that the 
leading logs L, especially the ones multiplied by a power of ln(M|/r|), are taken into account 
completely. As stated above, all sub-leading terms can not be generated neither by convolution 
nor by exponentiation. Therefore, we can estimate the error of the third order contribution to 
be about 1/L (actually we take 1/5) of its own size. For the third order hadronic contribution, 
derived by means of convolution, we should take into account also the intrinsic uncertainty in 
hadronic pairs, discussed above. As concerning the fourth and higher orders, we suggest to 
take the error just of the order O {a^L"^) e^e~ contribution (for all flavours together). 

The situation is better for the case of the final state pair corrections: the phenomenological 
treatment of hadronic pairs and the omitted effect of the interplay between the FS photons 
and pairs can not lead to more than 10 percent error in the FS pair RC, which is already small 
itself. Remind also the interference between the initial and final pair production amplitudes, 
discussed above. Collecting the sources of theoretical uncertainties, involved into our approach, 
we conclude: i) aX LEPl energies the error is about 0.4 permil on the wings and 0.2 permil on 
the top of the resonance; ii) at LEP2, when the radiative return to the Z-peak is forbidden, 
we have a reasonably good error of about 0.1 permil (for the realistic cut 2;min = 0.7225); in) 
when the radiative return is allowed, the size of the error can reach 1 permil. 

The approach to pair corrections described here is implemented into the semi-analytical 
code ZFITTER v.6.30 jB]. A treatment of higher order pairs (with pair-photon interplay) 
is available also in TOPAZO 14 , GEANTLE KORALW and some other codes. Some 

comparisons of the approaches were done in Ref. 0. 
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